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Abstract 

Dijet production has been studied in neutral current deep inelastic e~^p scattering 
for 470 < < 20000 GeV^ with the ZEUS detector at HERA using an inte- 
grated luminosity of 38.4 pb"^ Dijet differential cross sections are presented in a 
kinematic region where both theoretical and experimental uncertainties are small. 
Next-to-leading-order (NLO) QCD calculations describe the measured differential 
cross sections well. A QCD analysis of the measured dijet fraction as a function of 
allows both a precise determination of as{Mz) and a test of the energy-scale de- 
pendence of the strong coupling constant. A detailed analysis provides an improved 
estimate of the uncertainties of the NLO QCD cross sections arising from the parton 
distribution functions of the proton. The value of as{Mz), as determined from the 
QCD fit, is as{Mz) = 0.1166 ± 0.0019 (stat.)+J5:S {exp.)toZH (th.). 
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1 Introduction 



Dijet production in neutral current (NC) e+p deep inelastic scattering (DIS) provides a 
rich testing ground of the theory of the strong interactions of quarks and gluons, namely 
quantum chromo dynamics (QCD). At leading order (LO) in the strong coupling constant, 
2+1 jet production^ in NC DIS proceeds via the QCD-Compton {y*q qg with V = 
7, and boson-gluon fusion {V*g qq) processes. Thus, the differential cross section 
for dijet production is directly sensitive to a^, which is the fundamental parameter of the 
theory. Selecting a phase-space region where the perturbative QCD (pQCD) predictions 
are least affected by theoretical uncertainties provides a compelling test of QCD and 
permits a precise determination of the strong coupling constant. 

In this paper, measurements of the differential cross sections for dijet production in NC 
DIS are presented and compared with next-to-leading-order (NLO) pQCD predictions 
after correction for hadronisation effects. The phase-space region is restricted to high 
values of the virtuality, Q'^, of the exchanged boson, 470 < Q"^ < 20000 GeV^. In 
this region, the experimental uncertainties on the reconstruction of both the positron 
and the hadronic final state are smaller than at lower Q^. In addition, the theoretical 
uncertainties due to the modelling of the hadronic final state, to the parton distribution 
functions (PDFs) of the proton and to the higher-order contributions are minimised. 
The comparison of the differential cross sections with the theoretical predictions of the 
underlying hard processes provides a test of the NLO QCD description of dijet production. 
The analysis takes fully into account the correlation between the value of as{Mz) used in 
the determination of the proton PDFs and that in the calculation of the partonic cross 
sections. Furthermore, a detailed study of the theoretical uncertainties has been carried 
out which includes the statistical and correlated systematic uncertainties from each data 
set used in the determination of the proton PDFs. The QCD analysis yields a precise 
determination of as{Mz) and its energy-scale dependence. The twelve-fold increase in 
integrated luminosity, in combination with the improved experimental and theoretical 
methods used here, produces a significantly more accurate determination of as{Mz) with 
respect to the previous ZEUS measurement 

2 Theoretical framework 
2.1 Kinematics 

For a given e^p centre-of-mass energy, ^/s, the cross section for NC DIS, e+p — > e+ + X, 
depends on two independent kinematic variables, which are chosen to be and the 

^ Hereafter we denote the proton remnant by 
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Bjorken scaling variable, xbj, where = —q^ and xbj = Q^I^^P ■ and P (g) is the 
four-momentum of the incoming proton (exchanged virtual boson, V^*). 

For dijet production in NC DIS at LO, three additional variables, ^, and $, are needed 
in order to describe the kinematics of the two outgoing massless partons which form the 
two jets. The fraction ^ of the proton four-momentum carried by the incoming parton 
is defined by ^ = Xsj ■ (1 + Ml2/Q'^)i where M12 is the invariant mass of the two final- 
state partons. The variable Zp is defined for each outgoing parton i {i = 1, 2) with 
four-momentum pi by 



P-q 2 > E..i,2-Br'(l-<:os9i*')' 

where 9* is the parton scattering angle in the incoming parton- centre-of-mass system 
and Zp^i -|- Zp^2 = 1- Experimentally, Zp can be determined from the energy, -E^^*, and polar 
angle, of each of the two reconstructed jets in the laboratory frame. The angle $ 
is the azimuthal angle between the parton and lepton scattering planes in the l^*-parton 
centre-of-mass system. 



2.2 Cross-section calculation 

According to the QCD-improved parton model, as quantitatively expressed by the factori- 
sation theorem of QCD and perturbation theory, a NC DIS differential cross section^ 
da, can be written as 

da = ^ {C}) ■ rf(3-a(xP,a^(/iij),/i^,/i|) ■ (1 + ^had), (1) 

which can be interpreted as follows. The cross section has the form of a convolution of the 
partonic hard cross sections, dda, with the PDFs, fa^Xj^jp), of the colliding proton, with 
respect to the fraction x of the proton four-momentum carried by the incoming parton. 
The partonic cross sections describe the short- distance structure of the interaction and 
are calculable as power-series expansions in the strong coupling constant, which depends 
on the renormalisation scale, hr. The PDFs contain the description of the long-distance 
structure of the incoming proton. Their evolution with the factorisation scale, hf, at 
which they are probed follows the DGLAP equations 0. In Eq. (|I]), the dependence of 
the PDFs on the value of the strong coupling constant assumed in the DGLAP equations 
and on the parameters (collectively indicated by {(}) needed to model the x dependence 

^ The same symbol, dcr, is used for both the DIS inclusive, d(7tot, and the dijet, d(T2+ij differential cross 
sections. The same symbols with the superscript NLO indicate only the perturbative component (the 
first term on the right-hand side of Eq. (0)). 
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of the PDFs is explicitly indicated. In this analysis, the implicit as{Mz) dependence of 
the PDFs in the calculation of the NLO differential cross sections has been taken into 
account. The non-perturbative contribution, ^had, to the inclusive DIS cross section, atot, 
can be safely neglected in the high-Q^ region studied here. For the differential dijet cross 
sections, (5had was estimated using Monte Carlo (MC) models for the parton cascade and 
fragmentation (see Section 4). 

In the past few years there has been considerable theoretical progress on the development 
of general (i.e. process- and observable- independent) algorithms that allow a complete 
analytical cancellation of the soft and collinear singularities encountered in the calculation 
of NLO jet cross sections 0,0; flexible programs that compute arbitrary infrared- and 
coUinear-safe DIS observables in NLO QCD are now available 

3 Experimental setup 

The data sample was collected with the ZEUS detector ^ during the 1996-1997 data- 
taking period at HERA and corresponds to an integrated luminosity of 38.4 ± 0.6 pb~^. 
During this period HERA operated with protons of energy Ep = 820 GeV and positrons 
of energy = 27.5 GeV. 

The compensating uranium-scintillator calorimeter (CAL) covers 99.7% of the total 
solid angle. It is divided into three parts with respect to the polar angl^, 9, as viewed 
from the nominal interaction point: forward (FCAL, 2.6° < 6 < 36.7°), barrel (BCAL, 
36.7° < 9 < 129.1°), and rear (RCAL, 129.1° < 9 < 176.2°). Each section is sub- 
divided into towers which subtend solid angles between 0.006 and 0.04 steradian. Each 
tower is longitudinally segmented into an electromagnetic and one (RCAL) or two (FCAL, 
BCAL) hadronic sections. The electromagnetic section of each tower is further sub-divided 
transversely into two (RCAL) or four (BCAL, FCAL) cells. Under test-beam conditions, 
the calorimeter single-particle relative resolutions were 18%/ E (GeV) for electrons and 
35%/ a/ E (GeV) for hadrons. 

Tracking information is provided by the central tracking detector (CTD) |ll[ operating 
in a 1.43 T solenoidal magnetic field. The interaction vertex is measured with a typ- 
ical resolution along (transverse to) the beam direction of 0.4 (0.1) cm. The CTD is 
used to reconstruct the momenta of tracks in the polar angle region 15° < 9 < 164°. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as = — In (tan |), where the polar angle, 9, is measured with respect to the proton beam 
direction. 
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The transverse momentum, pt, resolution for full-length tracks can be parameterised as 
(T{pt)/pt = 0.0058 pt © 0.0065 © 0.0014/pi, with pt in GeV. 



The luminosity is measured using the Bethe-Heitler reaction e~^p e^'jp |]T^. The 
resulting small-angle energetic photons are measured by the luminosity monitor, a lead- 
scintillator calorimeter placed in the HERA tunnel at Z = — 107 m. 



4 Monte Carlo models 



NC DIS events including radiative effects were simulated using the HERACLES 4.5.2 [|13 



MC program with the DJANG06 2.4 IIJ] interface to the hadronisation programs. HER- 
ACLES includes corrections for initial- and final-state radiation, vertex and propagator 
terms, and two-boson exchange. The QCD cascade is simulated using the colour-dipole 
model |]15[ including the LO QCD diagrams as implemented^ in ARIADNE 4.08 [0 and, 
as a systematic check of the final results, with the MEPS model of LEPTO 6.5 [|l^. Both 
MC programs use the Lund string model [jl9| of JETSET 7.4 for the hadronisation. 
To estimate the uncertainty caused by the modelling of hadronisation, events were also 
generated with the HERWIG 5.9 [pll] program, in which the fragmentation is simulated 



according to a cluster model [22 



The ZEUS detector response was simulated with a program based on GEANT 3.13 p3 |. 
The generated events were passed through the simulated detector, subjected to the same 
trigger requirements as the data, and processed by the same reconstruction and offline 
programs. 

The simulations based on ARIADNE and LEPTO give a very good description of the 
shape and magnitude of the measured distributions of the inclusive event sample as a 
function of xbj and Q^. The shape of most of the measured distributions of the di- 
jet sample is also adequately described by these two programs. LEPTO gives a better 
description of the jet pseudorapidity distributions, while ARIADNE gives a better de- 
scription of the shape of the dijet cross section as a function of Q^. ARIADNE was used 
to calculate the acceptance corrections; both programs give very similar results. 

In addition, samples of events were generated without either exchange or electroweak 
radiative events, allowing the measured cross sections to be corrected for these effects, 
which at present are not included in the NLO QCD programs described below. 



A modified treatment of parton radiation at high Q 16 1 was included 
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5 NLO QCD calculations 



The NLO QCD calculations used in the analysis are based on the program DISENT 
0. The calculations make use of a generalised version ||^ of the subtraction method 



2J] and are performed in the massless MS renormalisation and factorisation schemes. 
The number of flavours was set to 5, the scales /i/j and were both set to Q, and 
c^sifJ'R) was evaluated using the two- loop formula |2^. Various sets of proton PDFs, 



hereafter generically referred to as MBFIT PDFs0, resulting from a recent analysis p6 



were used. This analysis provides the covariance matrix of the fitted PDF parameters and 
the derivatives as a function of x and of the PDFs with respect to these parameters. 
These are necessary to propagate the statistical and correlated systematic uncertainties 
of each data set used in the NLO DGLAP fit to the NLO QCD differential cross sections. 
The DISENT predictions have been cross-checked with results obtained with the program 
DISASTER-|--|- 0. Agreement at the 1 — 2% level has been found between these programs 
for all the observables discussed in this paper. 

The NLO pQCD predictions to be compared with the data were corrected by a bin-by-bin 
procedure for hadronisation effects according to da2+i = rfcr^*^ • C^^. The hadronisation 
correction factor, Chad (= (1 +^had)~^)5 "was defined as the ratio of the dijet cross sections 
before and after the hadronisation process, Chad = "cr2+i /"'^2+i ■ ^'^^ value 
of Chad was taken as the mean of the ratios obtained using the ARIADNE and HERWIG 
programs, for which the predictions agree typically within 5%. The hadronisation correc- 
tion was below 10% for most of the phase space. These programs only include the LO 
matrix elements for the partonic cross sections and higher-order contributions are simu- 
lated using various approximations to the parton cascade. This procedure for applying 
hadronisation corrections to the NLO QCD predictions was verified by checking that the 
shapes of the differential dijet cross sections at the parton level are reproduced by the 
NLO QCD calculations at the 10% level. 



6 Event selection and jet search 
6.1 The NC DIS data sample 

High-Q^ NC DIS events were selected according to criteria similar to those described in 
a recent ZEUS publication The events are characterised by a high-energy isolated 
positron in the detector, making them easy to distinguish from quasi-real photoproduction 
(Q^ ~ 0) and beam-gas interaction backgrounds. 

^ The MBFITIM set of proton PDFs assumes as{Mz) 0.118. 
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In addition to the energy and polar angle of the positron^ the variables 6, net transverse 
momentum, Pt, and total transverse energy, Et, were used for event selection, where 
6 = E — Pz, E is the total energy as measured in the CAL and P = Eirj. The sum 
runs over all calorimeter energy deposits Ei, and r,- is a unit vector along the line joining 
the reconstructed vertex and the geometric centre of cell i. The following criteria were 
applied offline: 

• to ensure that event quantities can be accurately determined, a reconstructed vertex 
with —50 < Z < 50 cm was required; 

• to suppress photoproduction events, in which the scattered positron escapes through 
the beam hole in the RCAL, 5 was required to be greater than 38 GeV. For perfect 
detector resolution, 6 is twice the positron beam energy (55 GeV) for fully contained 
DIS events, while, for photoproduction events, 6 peaks at much lower values. This cut 
also rejected events with hard initial-state QED radiation. The additional requirement 
(5 < 65 GeV removed cosmic ray background; 

• positrons were identified based on calorimeter cluster quantities and tracking: 

— to ensure high purity, the positron was required to have an energy of at least 
10 GeV; 

— to reduce background, isolated positrons were selected by requiring no more than 
5 GeV in calorimeter cells not associated with the scattered positron in an 77 — 
cone of radius 0.8 centred on the positron direction; 

— each positron with > 17.2° was required to match to a charged track of at least 
5 GeV momentum. For positrons beyond the tracking acceptance {6^ < 17.2°), 
the tracking requirement in the positron selection was replaced by a cut on the 
transverse momentum of the positron > 30 GeV and by the requirement 6 > 
44 GeV; 

— a fiducial volume cut was applied to the positron position which excludes the 
upper part of the central RCAL area (20 x 80 cm^) occluded by the cryogenic 
supply for the solenoid magnet. The transition regions between FCAL and BCAL, 
and BCAL and RCAL, corresponding to polar angles of 35.6° < 9 < 37.3° and 
128.2° < 9 < 140.2°, respectively, were also excluded; 

• to reduce further the background from photoproduction, the variable y = Q'^/{xBjs), 
estimated from the positron energy and angle, Ue, was required to satisfy i/e < 0.95; 

^ The positron polar angle, 9e, was determined from the associated track if the positron cluster was 
within the CTD acceptance, and otherwise from the position of the cluster and the reconstructed 
vertex. 
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• the net transverse momentum, Pj- , is expected to be close to zero and was measured 
with an error approximately proportional to ^^/EriGeV). To remove cosmic rays and 
beam-related background, PT(GeV) was required to be less than A^/ET{G(^^■ 

The kinematic variables Q^, xbj and y were determined using the double angle (DA) 



reconstruction method and are referred to as Ql,^, Xqa and y^A, respectively. The 
DA method is insensitive to errors in the absolute energy scale of the calorimeter and 
MC studies have shown it to be superior to other methods in the high-Q^ region 



considered here. 



6.2 Jet search and reconstruction 

The longitudinally invariant /c^-cluster algorithm was used in the inclusive mode to 



reconstruct jets in the hadronic final state both in data and in MC simulated events. 
At the detector level, the algorithm was applied to the energy deposits in the CAL cells 
after excluding those associated with the scattered positron candidate. The jet search 
was also applied to the generated hadrons and to the partons in the MC samples, as 
well as to the partons of the NLO QCD programs. The jet search was performed in the 
Breit frame |30|, in which the exchanged virtual boson is purely space-like with three- 



momentum q = (0, 0, —Q). In this frame, the selection of the jets in terms of transverse 
energy with respect to the direction of the virtual boson allows a natural suppression 
of the contribution due to single-jet events and strongly reduces the contamination from 
the proton remnant. The boost was performed using the four- momenta of the incoming 
and scattered positrons. The scattered positron four-momentum was calculated using the 
polar and azimuthal angles of the positron trackQ and the DA positron energy, defined as 

<DA = QDA/(2i5e(l + COS^e)). 

The jet algorithm uses quantities defined in the Breit frame and with respect to the 
direction of the incoming proton: the transverse energy, E^^, the pseudorapidity, tj'q, 
and the azimuthal angle, 0^, of the object i. The jet variables were defined according to 
the Snowmass convention |l3l| : 



7?jet _ \^ pi jet _ ^T,B Vb ,jet _ Zli ^T,B B 

^T,B - Z^^T,B 5 Vb - T^t ' ^B - • 

-^T,B ^T,B 



The comparison of the reconstructed jet variables between jets of hadrons and jets of 
CAL cells in MC generated events showed no significant systematic shift in the angular 



^ For events with 9e < 17.2°, the polar and azimuthal angles were calculated from the position of the 
positron calorimeter cluster and the reconstructed vertex position. 
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variables rf^ and (/>''J . However, the jet transverse energy as measured by the CAL under- 
estimated that of the jet of hadrons by an average of ~ 15%. This effect is due mainly to 
energy losses in the inactive material in front of the CAL and was corrected using 
the samples of MC generated events. 



6.3 Selected kinematic region 

The differential cross sections and dijet fraction presented in Section 7 are quoted for 
the DIS kinematic region defined by 470 < Q'^ < 20000 GeV^ and < y < 1, and 
for the jet selection criteria E)^^^ > 8 GeV, E}^^^ > 5 GeV and —1 < rf^lf^ < 2, 
where E^r^^^ (E^^^^) is the transverse energy of the jet in the Breit frame with the 
highest (second highest) transverse energy in the event. Only events with exactly two jets 
passing the above selection cuts were used to compute the cross sections. The final sample 
contains 1637 dijet events. The restriction to a high-Q^ region was chosen to avoid the 
large renormalisation-scale dependence of the NLO QCD dijet cross sections at lower 
(20 — 50% in the region 10 < < 100 GeV^) ^Sf, to reduce the uncertainty due to the 
proton PDFs (in particular the gluon density), and to improve the reconstruction of the 
boost to the Breit frame. The asymmetric cuts on the -E^*^ of the jets avoided infrared 
sensitive regions where the behaviour of the cross section as predicted by the NLO QCD 
programs is unphysical p4| . 



7 Results 



7.1 Dijet differential cross sections 

The cross sections measured in the kinematic region defined in the previous section were 
corrected for detector effects, QED radiative effects and Z°-exchange processes. In the 
presented dijet cross sections, the two jets were ordered according to decreasing pseudora- 
pidity in the Breit frame [rj^^'^ > Vs*'^)- The measurements^ of the differential dijet cross 
sections are presented in Figs. |l| and ^. The inclusive and dijet cross sections, as well as 
the dijet fraction, i?2+i(<5^) = {do-2+i/dQ'^)/ (datot/dQ^), are presented as a function of 
in Fig. |. 

In Figs. |T]to ^ the measured differential dijet cross sections and dijet fraction are compared 
to DISENT NLO QCD predictions, corrected for hadronisation effects. The hadronisation 



Tables of the results and their associated uncertainties are available in electronic form from the ZEUS 



WWW site at http://www-zeus.desy.de/zeus_papers/zeus-papers.htm] . They can also be obtained by 



contacting the authors. 
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correction, Chad, as well as its uncertainty, and the NLO QCD cross sections without 
hadronisation corrections are shown in Figs. |l| to ^. The QCD predictions, which assume 
as{Mz) = 0.118, provide a good overall description of both the shape and the magnitude 
of the measured cross sections. The only case where the agreement between data and 
NLO QCD is not good is the cross section as a function of the pseudorapidity of the most 
backward jet in the Breit frame (see Fig. |^d), where the data lie above the theoretical 
predictions in the region rf^^''^ > 0.5. However, hadronisation effects are particularly large 
for this variable and vary more rapidly than for other variables. The impact of such a 
discrepancy in the determination of as{Mz) is minor since the contribution to the total 
dijet cross section from the region f]^^^''^ > 0.5 is small. The dijet fraction increases with 
increasing Q"^ due to phase-space effects. For the cross sections as a function of the jet 
transverse energies and Q^, there is agreement at the ~ 10% level between data and 
theory over four orders of magnitude, demonstrating the validity of the description of the 
dynamics of dijet production by the NLO QCD hard processes. 



7.2 Experimental and theoretical uncertainties 



A detailed study of each of the main sources contributing to the systematic uncertainties 
of the measurements has been performed [^. These sources, for which a typical value of 
the systematic uncertainty in the dijet cross section da2+i/dQ'^ is indicated in parentheses, 
are listed below: 



• uncertainties in the positron identification efficiency and in the positron energy-scale 

m (1%); 

• uncertainties in the reconstruction of the boost to the Breit frame (1%); 

• use of the LEPTO program instead of ARIADNE to evaluate the acceptance correc- 
tions to the observed dijet distributions (2%); 

• uncertainties in the simulation of the trigger and in the cuts used to select the data 
(2%); 

• uncertainty of ±2% in the absolute jet energy scale of the CAL (±3.5%) obtained by 
examining the transverse momentum balance between the scattered positron and the 
jets. 



The last uncertainty is the dominant source of experimental systematic uncertainty and 
is strongly correlated between measurements at different points. It is shown as a light 
shaded band in Fig. In addition, there is an overall normalisation uncertainty of 1.6% 
from the luminosity determination, which is not included in the figures. 
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The NLO QCD predictions for the dijet cross sections are affected by the following: 

• uncertainties in the hadronisation correction, which were estimated as half the spread 
between the Chad values obtained using the ARIADNE (string) and HERWIG (cluster) 
models (1%); 

• uncertainties due to terms beyond NLO, which were estimated by varying /i^ between 
(5/2 and 2(5, keeping [Xp fixed at Q (6%); 

• uncertainty in the value of as(M^), which was estimated by repeating the calculations 
using MBFIT proton PDFs determined assuming a,(Mz) = 0.113 and 0.123 P§ (6%); 



uncertainties in the proton PDFs. The uncertainty due to the statistical and corre- 
lated systematic experimental uncertainties of each data set used in the determination 
of the MBFIT PDFs was calculated making use of the covariance matrix provided p6 



To estimate the uncertainties on the NLO QCD differential cross sections due to the 
theoretical uncertainties affecting the extraction of the NLO PDFs, the calculation 
of all the differential cross sections was repeated using a number of different MB- 



FIT PDFs obtained under different theoretical assumptions in the DGLAP fit p6 
This uncertainty was added in quadrature to that estimated above to give the total 
PDF-related uncertainty. The importance of taking the correlations among the PDF 
parameters into account is illustrated in Fig. ^, which shows, as a function of (5^, the 
relative uncertainties arising from the PDFs on the inclusive and dijet cross sections, 
and the dijet fraction. The uncertainties were obtained by including or ignoring the 
information provided by the covariance matrix mentioned above. The small uncer- 
tainty on i?2+i(Q^) due to the PDFs makes this observable particularly suitable for 
extracting (2.5% on datot/dQ'^, 4% on da2+i/dQ'^ and 1.5% on i?2+i((5^))- 



The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. 



7.3 Determination of a,^ 

The measured dijet fraction as a function of (5^, R2+i{Q'^)^ was used to determine as{Mz). 
The sensitivity of the measurements to the value of as{Mz) is illustrated in Fig. |^, which 
compares the measured -R2+i(<5^) with the NLO QCD calculations for three values of 
Ois{Mz). 

The procedure to determine as{Mz) was as follows: 
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NLO QCD calculations of i?2+i(Q^) were performed for three sets of MBFIT proton 
PDFs obtained assuming as{Mz) = 0.113, 0.118 and 0.123, respectively ||26|. The 



value of as{Mz) used in each partonic cross-section calculation was that associated 
with the corresponding set of PDFs; 

• for each bin, i, in Q^, the NLO QCD calculations mentioned above, corrected for 
hadronisation effects, were used to parameterise the as dependence of the dijet fraction 
according to the functional form: 

Rl_,,{as{Mz)) = A\ ■ as{Mz) + ■ al{Mz). (2) 

This parameterisation allows a simple description of the as dependence of i?2+i(Q^) 
over the entire range spanned by the MBFIT PDF sets, while using only three 
NLO calculations of the dijet fraction; 

• the value of as{Mz) was then determined by a x^-fit of Eq. (Q) to the measured 
-R2+i(<5^) values. 

This procedure correctly handles the complete a^-dependence of the NLO differential 
cross sections (the explicit dependence coming from the partonic cross sections and the 
implicit one coming from the PDFs) in the fit, while preserving the correlation between 
as and the PDFs. Its stability was checked with respect to variations in the PDFs and 
as-, as well as to alternative parameterisations of R2+i{Q'^ ■,as{Mz))■ 
^aking into account only the statistical errors on the measured dijet fraction, as{Mz) 
is found to be 0.1166 ± 0.0019 (stat.). The uncertainty on the value of as{Mz) due to 
the experimental systematic uncertainties of the measured dijet fraction was evaluated by 
repeating the analysis above for each systematic check. For the ratio R2+i{Q'^)i some of 
the systematic uncertainties largely cancel. The total experimental systematic uncertainty 
on the value of as{Mz) is t^ooss. The largest uncertainty is that on the jet energy scale. 

The following sources of theoretical uncertainties (evaluated as described in Section 7.2) 
and cross-checks on the extracted value of as{Mz) were considered ||33|| : 



terms beyond NLO: Aa,(M^) = 1°:°°^; 
uncertainties in the proton PDFs: Aa;s(M^) = t^oou; 
hadronisation effects: Aas(M^) = ±0.0005; 

the dependence on the use of a renormalisation scale that involves the jet variables. 
The analysis was repeated with = [E^^^^ + E^r^^j^) . The result is as{Mz) = 0.1125± 
0.0018 (stat.), which is almost equal to that obtained with fiR = Q/2. This source of 
uncertainty was therefore neglected; 
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• the fit procedure. This was cross-checked by repeating the as{Mz) determination 
using the five sets of proton PDFs of the CTEQ4 "A-series" ^5] and the three MRST 
sets, central, as]] and ii The results are in good agreement with the central 
value determined above. This source of uncertainty was therefore neglected 



The total theoretical uncertainty of ^0^0044 was obtained by adding in quadrature the 
uncorrelated uncertainties on as{Mz) due to the first three items mentioned above. 

The value of as{Mz) as determined from the measured -R2+i(Q^) is therefore: 

as{Mz) = 0.1166 ±0.0019 (stat.)+S:S (exp.)+S:E (th.). 
This result is consistent with the current PDG world average, as{Mz) = 0.1181 ± 0.0020 



2^], a review from Bethke p^, and recent determinations by the HI Collaboration [R8 



7.4 The energy scale dependence of as 

A consistency test for the scale dependence of the renormalised strong coupling constant 
predicted by the renormalisation group equation was carried out by repeating the QCD fit 
of the dijet fraction in five bins. The principle of the fit is the same as outlined above, 
with the only difference being that the dependence of the dijet fraction in Eq. 
was parameterised not in terms of as{Mz) but in terms of as{{Q)), where (Q) is the 
mean value of Q in each bin. The measured as{{Q)) values, with their experimental and 
theoretical systematic uncertainties estimated as for as{Mz), are shown in Fig. ||. The 
measurements are compared with the renormalisation group predictions obtained from 
the PDG as{Mz) value and its associated uncertainty. The values are in good agreement 
with the predicted running of the strong coupling constant over a large range in Q. 



8 Summary 

Differential dijet cross sections have been measured in neutral current deep inelastic e~^p 
scattering for 470 < < 20000 GeV^ with the ZEUS detector at HERA. These mea- 
surements, which use events with only two jets, were performed in a kinematic region 
where both theoretical and experimental uncertainties are small. Next-to-leading-order 
QCD calculations give a good description of the shape and magnitude of the kinematic 
distributions and cross sections. For the cross sections as a function of the jet transverse 
energies and Q^, there is agreement at the ~ 10% level between data and theory over four 
orders of magnitude, demonstrating the vahdity of the description of the dynamics of dijet 
production by the NLO QCD hard processes. A QCD fit of the measured dijet fraction 
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as a function of provides both a precise determination of the strong couphng constant 
and a test of its energy-scale dependence. A comprehensive analysis of the uncertainties 
of the calculations has been carried out, which takes into account the dependence of the 
proton PDFs on the assumed value of as and the statistical and correlated systematic 
uncertainties from each data set used in the determination of the proton PDFs. The 
value of as{Mz) as determined by fitting the next-to- leading-order QCD calculations to 
the measured dijet fraction is 

as{Mz) = 0.1166 ±0.0019 (stat.)+;;S (exp.)^S:E (^h-) , 



in good agreement with the world average E^. The value of function of Q is in 



good agreement, over a wide range of Q, with the running of as predicted by QCD. 
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Zp^i range 


da2+i/dzp^i [pb] 


Astat 


Asyst 


Aes 


0.00 - 0.05 


1.118 ■ 10^ 


(±0.244) ■ 10^ 

V / 


(^0 09?) ■ 
v— u.uy 1 / 


/+0.089N inl 
V — U.Uoy / 


0.05 - 0.15 


9.777 ■ 10^ 


(±0.555) • 10^ 

V / 


(+0 5?9) ■ 


(^oioo) ■ 10' 


0.15 - 0.25 


1.233 ■ 10^ 


(±0.063) ■ 10^ 

V / 


(+0 068) ■ 


("■ooi) ■ 10' 


0.25 - 0.35 


8.469 ■ 10^ 


(±0.520) • 10^ 


(^0 482) ■ 10' 


(^0 ?58) • 10' 


0.35 - 0.45 


7.738 ■ 10^ 


(±0.505) • 10^ 


/+0.425N . 


doil") ■ 10' 


0.45 - 0.55 


5.473 ■ 10^ 


(±0.438) ■ 10^ 


/+0.288N inl 
\-0.U5) ' 


/+0.125N . iQl 
U0.I33J 


0.55 - 0.65 


5.259 ■ 10^ 


(±0.453) ■ 10^ 


it'oll) ■ 10' 


/+0.138N . iQl 
1-0.148J 


0.65 - 0.75 


3.938 ■ 10^ 


(±0.415) ■ 10^ 


/+0.240N . .Ql 
I-O.451J 


/+0.167N . .Ql 
I-O.I35J 


0.75 - 0.85 


2.064 ■ 10^ 


(±0.349) ■ 10^ 


/+0.076N . .Ql 
1-0.416J 


/+0.135N . -|Ql 
I-O.I37J 


0.85 - 0.95 


4.436 


±1.537 


+1.261 
-1.232 


+0.522 
-0.325 



Table 1: The differential dijet cross section da2+i/ dzp^i. For each bin in Zp^i, 
the measured cross section, the statistical uncertainty, Agtat, and the systematic 
uncertainty (not) associated with the energy scale of the jets, Aes (Asyst), are given. 
The overall normalisation uncertainty of 1.6% due to the luminosity determination 
is not included. 



logiol^^Bj) range 


da2+i/d\og^Q{xBj) [pb] 


Astat 


Agyst 


Aes 


-2.20 - -2.00 


2.681 ■ 10^ 


(±0.267) ■ 10^ 


/+0.24U . .Ql 
I-O.O93J 


■ 10' 


-2.00 - -1.80 


4.794 ■ 10^ 


(±0.286) ■ 10^ 


/+0.139N . .Ql 
1-0.118J 


(;S:}}^) ■ 10' 


-1.80 - -1.60 


5.902 ■ 10^ 


(±0.312) ■ 10^ 


/+0.143^ . .Ql 
I-O.472J 


/+0.162N . -|Ql 

I-0.166J 


-1.60 - -1.40 


5.821 • 10^ 


(±0.319) ■ 10^ 


/+0.110N . .pi 

I-O.497J -"-U 


/+0.175N . -|Ql 

I-0.180J 


-1.40 - -1.20 


4.105 ■ 10^ 


(±0.263) ■ 10^ 


/+0.341N . iQl 
I-O.I41J 


/+0.146N . iQl 
U0.138J 


-1.20 - -1.00 


2.750 ■ 10^ 


(±0.203) ■ 10^ 


n'^) ■ 10' 


/+0.095N inl 
1-0.070^ ■ 


-1.00 - -0.80 


1.195 ■ 10^ 


(±0.131) ■ 10^ 


(1^:1?^) ■ 10' 


/+0.04U . .Ql 
I-O.O34J 


-0.80 - -0.60 


5.078 


±0.793 


+0.323 
-0.670 


+0.176 
-0.105 


-0.60 - -0.40 


1.373 


±0.414 


+0.363 
-0.552 


+0.042 
-0.023 



Table 2: The differential dijet cross section da2+i/ dlogiQ^XBj)- Other details 
are as described in the caption to Table |1|. 
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logiolO range 


da2+i/dlogio(0 [pbj 


A 


A 

^syst 




-2.1875 - -1.8750 


3.386 


±0.939 


+U.y ( U 
-0.736 


-0.078 


-1.8750 - -1.5625 


2.726 ■ 10^ 


(±0.179) ■ 10^ 


1-0.045^ 


/+0.021N , 1 Ql 


-1.5625 - -1.2500 


6.379 ■ 10^ 


(±0.266) ■ 10^ 


/_Ln nfici \ ^ ^ 1 
1-0.334^ ■ 


/_i_n 1 /I v \ ^1 
1-0.136J 


-1.2500 - -0.9375 


5.997 ■ 10^ 


(±0.260) ■ 10^ 


/+0.156X . iQl 
1-0.224/' -"-'J 


/+0.199N inl 
I-O.I89J ■ 


-0.9375 - -0.6250 


2.232 ■ 10^ 


(±0.142) ■ 10^ 


/+0.105N . iQl 
I-O.I54J 


/+0.126N . iQl 


-0.6250 - -0.3125 


2.434 


±0.400 


+0.305 
-0.229 


+0.177 
-0.123 



Table 3: The differential dijet cross section (icr2+i/'^ logio(0- Other details are 
as described in the caption to Table |^. 



Mjj range [GeV] 


da2+i/dMjj [pb/GeV] 


Astat 


Asyst 


Aes 


10.00 - 16.00 


7.462 ■ 10-1 


(±0.642) ■ 10-1 


/+0.562N ir|-l 
U1.434J ■ -"-"J 


/+0.016N . 1Q-I 
I-O.OO5J 


16.00 - 22.00 


2.404 


±0.122 


+0.055 
-0.055 


+0.034 
-0.035 


22.00 - 28.00 


2.127 


±0.109 


+0.048 
-0.041 


+0.057 
-0.053 


28.00 - 35.00 


1.272 


±0.077 


+0.060 
-0.038 


+0.046 
-0.042 


35.00 -45.00 


7.494 ■ 10-1 


(±0.497) ■ 10-1 


/+0.090N -|n-l 
I-O.309J ■ -"-U 


/+0.256N in-i 
1-0.246J ■ 


45.00 - 60.00 


3.430 ■ 10-1 


(±0.284) ■ 10-1 


/+0.100N . 1Q-I 
1-0.484^ 


/+0.270N . 1Q-I 
1-0.195^ 


60.00 - 80.00 


9.605 ■ 10-2 


(±1.272) ■ 10-2 


/+0.326N 10-2 
I-O.419J ■ 


/+0.619N in-2 
I-O.709J ■ 


80.00 - 120.00 


1.578 ■ 10-2 


(±0.353) ■ 10-2 


/+0.122N . -|Q-2 
I-O.219J 


/+0.117N . -|Q-2 
I-O.I45J 



Table 4: The differential dijet cross section da2+i/dMjj. Other details are as 
described in the caption to Table |l|. 



E^*^^ range [GeV] 


da^+i/dE':^"^ [pb/GeV] 


Astat 


Asyst 


Aes 


5.00 - 8.00 


2.956 


±0.175 


+0.154 
-0.243 


+0.072 
-0.081 


8.00 - 12.00 


5.977 


±0.238 


+0.105 
-0.201 


+0.094 
-0.090 


12.00 - 16.00 


2.736 


±0.154 


+0.039 
-0.033 


+0.083 
-0.081 


16.00 - 20.00 


1.327 


±0.103 


+0.055 
-0.030 


+0.072 
-0.026 


20.00 - 24.00 


8.165 ■ 10-1 


(±0.801) • 10-1 


/+0.096N in-l 
1-0.397^ ■ 


/+0.245N 1Q-I 
1-0.566J 


24.00 - 32.00 


2.848 ■ 10-1 


(±0.333) ■ 10-1 


{t'oTi) ■ 10-' 


/+0.236N -in-i 

I-0.183J ■ 


32.00 - 50.00 


5.960 ■ 10-2 


(±0.879) • 10-2 


/+0.997N in-2 
1-0.449^ ■ 


/+0.572N 1Q-2 
I-O.497J 


50.00 - 80.00 


3.991 • 10-3 


(±1.203) ■ 10-3 


/+1.086N in-3 
I-O.902J ■ 


/+0.548N in-S 
1-0.472^ ■ 



Table 5: The differential dijet cross section da2+i/dEi^^ . Other details are as 
described in the caption to Table 
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jet.l 

Vb range 


da 2+1 /drj^^^'^ [pb] 


A 


A 

^syst 


A 

^ES 


-1.00 - -0.50 


3.051 


±0.419 


+0.360 
-0.277 


+0.155 
-0.141 


-0.50 - 0.00 


1.170 • 10^ 


(±0.084) • 10^ 


i-o ill) ■ 10' 


/-LO 060\ ^ ^1 

U0.051J ' 


0.00 - 0.50 


3.260 ■ 10^ 


(±0.154) ■ 10^ 


(^OAOs) ■ 10^ 


(-aog?) ■ 10^ 


0.50 - 1.00 


3.204 ■ 10 


(±0.151) ■ 10^ 




/J-O 068\ 1 n1 

I-O.O75J ■ '-^ 


1.00 - 1.50 


2.054 ■ 10^ 


(±0.119) ■ 10^ 


i'+0.083n ml 
U0.IO4J ■ 


i'+0.046n 1 nl 
1-0.049/ ■ 


1.50 - 2.00 


9.035 


±0.748 


+0.570 
-0.742 


+0.229 
-0.194 


2.00 - 2.50 


2.561 


±0.410 


+0.491 
-0.324 


+0.061 
-0.058 



Table 6: The differential dijet cross section da 2+1 / drf^*"'^ . Other details are as 
described in the caption to Table [I|. 



Ei^*'^ range [GeV] 


da2+i/dE^T^f [pb/GeV] 


Astat 


Asyst 


Aes 


5.00 - 8.00 


2.775 


±0.168 


+0.098 
-0.235 


+0.074 
-0.082 


8.00 - 12.00 


5.971 


±0.240 


+0.142 
-0.154 


+0.103 
-0.092 


12.00 - 16.00 


2.762 


±0.152 


+0.074 
-0.071 


+0.088 
-0.076 


16.00 - 20.00 


1.403 


±0.105 


+0.011 
-0.034 


+0.062 
-0.050 


20.00 - 24.00 


7.437 ■ 10-1 


(±0.763) • 10-1 


t'lir.) ■ 10-' 


('+0.397N . 1 Q-1 
1-0.406^ 


24.00 - 32.00 


2.995 ■ 10-1 


(±0.339) ■ 10-1 


n^r) ■ 10-' 


/+0.182N . -|Q-1 
1-0.198J 


32.00 - 50.00 


6.834 ■ 10-2 


(±0.966) ■ 10-2 


/+0.64U . -|Q-2 
1-0.239^ 


/+0.551N -|n-2 
I-O.452J ■ ^'^ 


50.00 - 80.00 


3.189 ■ 10-=^ 


(±1.205) ■ 10-3 


^+0.98^ m-3 
1-0.621/' ■ 


/+0.203N m-3 
1-0.251^ ■ 



Table 7: The differential dijet cross section da2+i/dE^^ . Other details are as 
described in the caption to Table |l|. 



jet,2 

7]-'^ range 


da2+i/dr]^j^^''^ [pb] 


Astat 


Agyst 


Aes 


-2.00 - -1.50 


2.507 


±0.458 


+0.219 
-0.670 


+0.202 

-0.190 


-1.50 - -1.00 


1.054 ■ lOi 


(±0.084) ■ lOi 


/+0.043N . .Ql 
1-0.084J 


/+0.069N ml 
1-0.067^ ■ 


-1.00 - -0.50 


3.328 ■ lOi 


(±0.151) ■ lOi 


/+0.080N inl 
U0.I70J ■ 


/+0.154N . Ql 
U0.I34J 


-0.50 - 0.00 


3.735 ■ lOi 


(±0.163) ■ lOi 


/+0.067>| ml 
U0.I30J ■ 


/+0.073N . 1 Ql 
1-0.066/ 


0.00 - 0.50 


1.885 ■ lOi 


(±0.112) ■ lOi 


/+0.062>, ml 
1-0.066J ■ 


/+0.012N . -|Ql 
I-O.OI4J 


0.50 - 1.00 


6.868 


±0.627 


+0.504 
-0.252 


+0.070 
-0.080 


1.00 - 1.50 


1.914 


±0.349 


+0.288 
-0.273 


+0.066 
-0.057 



Table 8: The differential dijet cross section da2+i/ drf^^'"^ . Other details are as 
described in the caption to Table |l|. 
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range [GeV^] 


da,otldQ^ [pb/GeV2] 


Astat 


Asyst 


470. - 800. 


1.217 


±0.011 


+0.015 
-0.006 


800. - 1500. 


2.800 ■ 10-1 


(±0.033) ■ 10-1 


/+0.023N in-l 

l-o.oigJ ■ -"-u 


1500. - 2500. 


6.509 ■ 10-2 


(±0.130) • 10-2 


/+0.103N -|r|-2 
l-0.090J ■ 


2500. - 5000. 


1.230 ■ 10-2 


(±0.036) • 10-2 


/+0.022\ -|r|-2 
I-O.O20J ■ 


5000. - 20000. 


7.378 • 10-4 


(±0.400) • 10-4 


/+0.051N in-4 
1-0.163/' ■ 



Table 9: The differential inclusive cross section da tot I dQ'^. Other details are as 
described in the caption to Table |l|. 



g2 range [GeV2] 


da2+i/dQ^ [pb/GeV2] 


Astat 


Agyst 


Aes 


470. - 800. 


8.061 ■ 10-2 


(±0.306) ■ 10-2 


/+0.212N . -|Q-2 
I-O.294J 


/+0.287A . 1Q-2 
1-0.256J 


800. - 1500. 


2.415 ■ 10-2 


(±0.105) ■ 10-2 


/+0.039N in-2 
I-O.O61J ■ 


I-O.O72J -"-U 


1500. - 2500. 


6.802 ■ 10-3 


(±0.453) ■ 10-3 


/+0.300N in-3 

1-0.343J ■ 


/+0.159N in-3 
I-O.I60J ■ 


2500. - 5000. 


1.637 ■ 10-=^ 


(±0.139) ■ 10-3 


/+0.046N ir,-3 
I-O.I24J ■ 


/+0.033N 1 n-3 
1-0.028^ ■ 


5000. - 20000. 


1.246 ■ 10-4 


(±0.170) ■ 10-4 


/+0.006>, in-4 

1-0.139J ■ 


/+0.013N in-4 

l-o.oisJ ■ 



Table 10: The differential dijet cross section do2+\l dQ^ . Other details are as 
described in the caption to Table |l|. 



g2 range [GeV2] 


-R2+1 


Astat 


Agyst 


Aes 


470. - 800. 


6.625 ■ 10-2 


(±0.245) ■ 10-2 


/+0.168\ in-2 
1-0.278^ ■ 


/+0.236N in-2 
I-O.2I0J ■ 


800. - 1500. 


8.627 • 10-2 


(±0.362) ■ 10-2 


/+0.144N . .Q-2 
1-0.267J 


/+0.255N in-2 
1-0.256J ■ 


1500. - 2500. 


1.045 ■ 10-1 


(±0.066) ■ 10-1 


/+0.043N . .Q-1 
I-O.O57J -"-U 


/+0.024N . .Q-l 
I-O.O25J 


2500. - 5000. 


1.331 ■ 10-1 


(±0.106) ■ 10-1 


/+0.039N in-l 

I-0.119J ■ -"-U 


/+0.027N . 1Q-I 
I-O.O23J -"-U 


5000. - 20000. 


1.688 ■ 10-1 


(±0.211) • 10-1 


/+0.008N in-i 
1-0.185J ■ 


/+0.017N . 1Q-I 
I-O.O24J 



Table 11: The dijet fraction i?2+i(<5^)- Other details are as described in the 
caption to Table 0. 
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Table 12: The QCD predictions for the differential dijet cross section as a function of Zp^i. For each bin in Zp^i the 
following quantities are given: the weighted mean value < Zp^i >, the pure NLO QCD cross section, the uncertainty 
due to the renormalization scale, A^^, the uncertainty due to ^^(M^), A^^, the uncertainty due to the proton PDFs 
(experimental), A^^J, the uncertainty due to the proton PDFs (theoretical), A*^|°, and the hadronisation correction, Chad, 
with its associated uncertainty, AChad- 



^ogioix Bj) range 


< logio(a;Bj) > 


d(^2+i^/d^og^oi^Bj) [pb] 








A thco 
^Pdf 


Chad ± AChad 


-2.2 - -2.0 


-2.08 


2.509 ■ 10^ 






(±0.140) ■ 10^ 


/+0.036\ inl 
1-0.040/' ■ ^'^ 


1.04± 0.02 


-2.0 - -1.8 


-1.90 


5.025 ■ 10^ 


/+0.285N ir|l 
1-0.297^ ■ 


mr^) ■ 10^ 


(±0.248) ■ 10^ 


/+0.064\ -|r,l 

1-0.072/' ■ 


1.06± 0.01 
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-1.70 


6.477 ■ 10^ 


/+0.336\ inl 
l-0.36lJ ■ 




(±0.256) ■ 10^ 


/+0.076\ inl 
1-0.078/' ■ 


1.08± 0.02 


-1.6 - -1.4 


-1.50 


6.403 ■ 10^ 


l-0.34lJ -LU 




(±0.205) ■ 10^ 


/+0.079\ 1 nl 
1-0.068/' ■ 


1.10± 0.01 


-1.4 - -1.2 
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1-0.083^ 


(±0.064) ■ 10^ 


/+0.017>, inl 
1-0.017/) -"-U 
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-0.8 - -0.6 


-0.72 


5.604 


+0.283 
-0.282 


+0.322 
-0.263 


±0.244 


+0.041 
-0.069 


1.11± 0.01 


-0.6 - -0.4 


-0.52 


1.650 


+0.088 
-0.082 


+0.061 
-0.051 


±0.062 


+0.000 
-0.023 


1.14± 0.05 



Table 13: The QCD predictions for the differential dijet cross section as a function of logiQ{xBj)- Other details are as 
described in the caption to Table 12. 



logio(0 range 


< logio(0 > 


d<r/c?logio(e) [pb] 


A^. 


Aa. 




A thco 
^pdf 


Chad ± AChad 


-2.1875 - -1.8750 


-1.94 


2.825 


+0.737 

-0.503 


+0.233 
-0.258 


±0.166 


+0.047 
-0.051 


1.18± 0.05 


-1.8750 - -1.5625 


-1.68 


2.951 • 10^ 


^+0.363\ inl 
1-0.290/' ■ 


/+0.167N . inl 
1-0.205/' 


(±0.180) • 10^ 


/+0.043\ 1 nl 
1-0.048/' ■ 


1.09± 0.01 
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-1.40 


6.736 • 10^ 


ailt) ■ 10^ 
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1-0.434/' ^'^ 


(±0.318) • 10^ 


/+0.075N . inl 
1-0.076/' ^'^ 


1.06± 0.01 


-1.2500 - -0.9375 


-1.11 


6.746 • 10^ 




/+0.453\ 1 nl 
1-0.448/' ■ 


(±0.273) • 10^ 
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1-0.069/' ■ 


1.11± 0.01 


-0.9375 - -0.6250 


-0.82 


2.569 • 10^ 


/+0.002\ inl 

1-0.059/' ■ 


it'oiTr) ■ 10' 


(±0.162) • 10^ 


/+0.041\ inl 

1-0.040/' ■ 


1.14± 0.02 


-0.6250 - -0.3125 


-0.53 


3.378 


+0.000 
-0.059 


+0.177 
-0.122 


±0.220 


+0.027 
-0.067 


1.11± 0.06 



Table 14: The QCD predictions for the differential dijet cross section as a function of logio(0- Other details are as 
described in the caption to Table 12. 
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Table 15: The QCD predictions for the differential dijet cross section as a function of Mjj. Other details are as 
described in the caption to Table 12. 

oo 
oo 



E^^^^ range [GeV] 


^ ^T,B ^ 


da^^\o/dE^:r'f [pb/GeV] 








A theo 
^pdf 


Chad ± AChad 


5. - 8. 


7.06 


3.722 


+1.174 
-0.790 


+0.429 
-0.397 


±0.108 


+0.056 
-0.062 


1.10± 0.01 


8. - 12. 


9.72 


6.637 


+0.023 
-0.176 


+0.304 
-0.366 


±0.407 


+0.075 
-0.065 


1.16± 0.03 


12. - 16. 


13.77 


2.905 


+0.000 
-0.067 


+0.145 
-0.163 


±0.190 


+0.034 
-0.031 


1.04± 0.00 


16. - 20. 


17.83 


1.420 


+0.000 
-0.030 


+0.079 
-0.083 


±0.096 


+0.018 
-0.019 


1.01± 0.01 


20. - 24. 


21.92 


7.071 ■ 10-1 


/+0.000N 1A-1 
1-0.218J ■ 


n^.) ■ 10-^ 


(±0.535) ■ 10-1 


/+0.105N . .A-l 
U0.IO7J -"-U 


1.00± 0.03 


24. - 32. 


27.28 


2.951 ■ 10-1 


/+0.001N . |A-1 

Uo.iisJ 




(±0.214) ■ 10-1 


(Toll) ■ 10-^ 


1.01± 0.03 


32. - 50. 


37.77 


5.730 ■ 10-2 


/+0.098N -|n-2 
l-0.38lJ ■ 


/+0.464N -|r|-2 
I-O.339J ■ 


(±0.482) ■ 10-2 




1.01± 0.03 


50. - 80. 


57.65 


3.908 • 10-3 


/+0.280N 1A-3 
1-0.609^ ■ 


/+0.324N 1A-3 
I-O.I62J ■ 


(±0.441) • 10-3 


/+0.094\ 1 A-3 
1-0.244^ ■ 


1.04± 0.03 



Table 16: The QCD predictions for the differential dijet cross section as a function of E^^^. Other details are as 
described in the caption to Table 12. 



rfg ' range 


. jet.l , 

<Vb > 


J NLO 1 J jet.l r 11 


A 


A 




A thfio 


(^had ± AOhad 


-1.0- -0.5 


-0.68 


5.156 


-j-U-oUO 
-0.294 


"hU-ooo 
-0.316 


±0.257 


-t-U.U04 

-0.046 


1.47± 0.05 


-0.5 - 0.0 


-0.21 


1.972 • 10^ 


/H-U.UDD\ -1 rvl 
V -0.086^ ■ 


/+0.115\ irvl 
V-O.I20J 


(±0.109) • 10^ 


/ _i_n no /I \ „ 1 
/+0.U24\ 1 rvl 

V-O.O22J 


1.34± 0.05 


0.0 - 0.5 


0.26 


3.596 • 10^ 


/+0.126\ 1 rvl 
V-0.163J 


(^0:222) • 10' 


(±0.197) • 10^ 


/+U.U40\ -1 rvl 
1-0.044^ ■ 


1.11± 0.01 


0.5 - 1.0 


0.73 


3.136 • 10^ 


V-0.168J 


/_l_n 1 Qyl \ - ^ 1 
V-O.2O5J -"-^ 


(±0.171) • 10^ 


/J-fl C\A 1 \ - ^ 1 

daoli) • 10' 


1.02± 0.02 


1.0-1.5 


1.22 


1.888 • 10^ 




/+0.129N . 
U0.I34J 


(±0.099) • W 


/+0.025\ -irvl 
1-0.026/* ■ 


0.99± 0.04 


1.5 - 2.0 


1.71 


8.064 


+0.884 
-0.724 


+0.596 
-0.617 


±0.413 


+0.106 
-0.112 


1.00± 0.01 


2.0 - 2.5 


2.18 


2.406 


+0.384 
-0.287 


+0.197 
-0.204 


±0.120 


+0.031 
-0.033 


1.08± 0.01 



Table 17: The QCD predictions for the differential dijet cross section as a function of 77^^*''. Other details are as 
described in the caption to Table 12. 



CO 



E^^*/ range [GeV] 


< E{^^ > 


da^^\^/dEiy [pb/GeV] 






Pdf 


^theo 
pdf 


C'had i AChad 


5. - 8. 


7.10 


3.623 


+1.148 
-0.771 


+0.425 
-0.390 


±0.109 


+0.056 
-0.063 


1.20± 0.02 


8. - 12. 


9.69 


6.737 


+0.053 
-0.196 


+0.315 
-0.376 


±0.398 


+0.076 
-0.066 


1.13± 0.02 


12. - 16. 


13.75 


2.877 


+0.000 
-0.062 


+0.139 
-0.159 


±0.185 


+0.033 
-0.031 


1.04± 0.01 


16. - 20. 


17.70 


1.367 


+0.000 
-0.033 


+0.071 
-0.076 


±0.094 


+0.017 
-0.016 


1.02± 0.02 


20. - 24. 


21.70 


7.527 • 10-1 


/+0.000\ -1 rv-1 
U0.I54J ■ 


/+0.464\ in-l 
{-0.455) ' 


(±0.527) ■ 10-1 


/+0.106N . 1Q-I 
1-0.108/* 


1.00± 0.03 


24. - 32. 


27.18 


3.020 • 10-1 


/+0.000\ -1 rv-1 
1-0.098^ ■ 


1-0.184J -"-U 


(±0.212) ■ 10-1 


/+0.052\ in-l 
I-O.O57J ■ -"-U 


0.99± 0.03 


32. - 50. 


37.71 


5.643 • 10-2 


/+0.117^, . -|Q-2 
I-O.4I0J 


/+0.446N ir,-2 
I-O.327J ■ 


(±0.476) • 10-2 


ariii) ■ 10-^ 


1.02± 0.03 


50. - 80. 


57.61 


3.867 • 10-=^ 


/+0.288>, -|r,-3 
U0.6I8J ■ 


/+0.315X . -|Q-3 
U0.I61J -"-U 


(±0.440) ■ 10-3 


(i°i) . 10-3 


0.97± 0.03 



Table 18: The QCD predictions for the differential dijet cross section as a function of E^^"^. Other details are as 
described in the caption to Table 12. 



jet,2 „ „^ 

r]^ range 


<Vb > 


J NLO 1 J j^t,2 r -i 1 


A 






A theo 


/~i _|_ A /-y 

(^had ± ^(^had 


-2.0- -1.5 


-1.69 


4.973 


-hU.o04 
-0.317 


_i_n Q/1 Q 
-0.309 


±0.264 


_i n ncyi 
-0.051 


1.51± 0.05 


-1.5 - -1.0 


-1.20 


1.972 • 10^ 


/+0.151\ lAl 
V-O.I35J -"-^ 


/+U.135\ 1 nl 
V-O.I33J 


(±0.104) • 10^ 


/ 1 n no /I \ _ ^ 1 
V-0.022y' ■ 


1.40± 0.00 


-1.0 - -0.5 


-0.74 


4.169 • 10^ 


/ -un 1 \ ^ 1 

(toirn) • 10' 


(lo:258) • 10' 


(±0.238) • 10^ 


/_i_n iT^i \ -i ^1 

(loiols) • 10' 


1.21± 0.03 


-0.5 - 0.0 


-0.28 


3.427 • 10^ 


/+0.132X , -,r,l 
V-0.164J -"-^ 


/+0.195X . -.rvl 
V-0.216J -"-^ 


(±0.180) • 10^ 


/+0.044\ -. r,l 
V-O.O43J ■ 


0.99± 0.02 


0.0-0.5 


0.20 


1.558 • 10^ 


/+0.100\ . -,r,l 


/+0.109X . -, Ql 

v-o.iioJ 


(±0.077) • 10^ 


/+0.023N 1 r,l 

V-0.024J ■ 


0.92± 0.01 


0.5 - 1.0 


0.69 


4.185 


+0.459 
-0.377 


+0.371 
-0.342 


±0.209 


+0.071 
-0.079 


0.96± 0.00 


1.0-1.5 


1.19 


9.857 • 10-^ 


■ 10-' 


/+1.041X . 
U0.945J 


(±0.454) • 10"^ 


/+0.170N . -iQ-l 
U0.201J 


1.36± 0.03 



Table 19: The QCD predictions for the differential dijet cross section as a function of 77^^*'^. Other details are as 
described in the caption to Table 12. 



Q2 range [GeV^] 




daf^^'IdC^ [plVGeV^] 






A Cxp 


.\ thoo 


170. 800. 


599. 


1.180 


+ 0.008 
-0.011 


+0,002 
-0.000 


±0.032 


+0.017 
-0.011 


800. - 1500. 


1059. 


2.756 • 10-^ 


/+0.016N . -|Q-1 
U0.O21J 


/+0.006\ -1 rv-l 

U0.012J ■ 


(±0.067) • 10-1 


/+0.034>, , .,^-1 
V-0.026J 


1500. - 2500. 


1890. 


6.282 • 10-2 


/+0.028\ -1 rv-2 
V-0.036/' ■ 


/+0.014\ 1 r,-2 
V-0.015^ ■ 


(±0.146) • 10-2 


(+0m4\ -ir,-2 
1-0.046^ ■ 


2500. - 5000. 


3363. 


1.257 • 10-2 


/+0.004\ -1 r,-2 
V-O.OO5J ■ -"-^ 


/+0.001\ -1 A-2 

v-o.oooJ ■ 


(±0.030) • 10-2 


/+0.009\ -irv-2 
V-O.OO7J 


5000. - 20000. 


7801. 


7.003 • 10-^ 


/+0.010\ -1 r,-4 
V-O.OI2J ■ 


/+0.057\ 1 r,-4 
V -0.035 J ■ 


(±0.188) • 10-^ 


/+0.003\ -ir,-4 
U0.O31J ■ 



Table 20: The QCD predictions for the differential inclusive cross section as a function of Q"^. Other details are as 
described in the caption to Table 12. 



range [GeV^] 




da^_^f/dQ^ [pb/GeV2] 






A cxp 


A thco 
pdf 


Chad ± AChad 


470. - 800. 


604. 


8.985 ■ 10-2 


/+0.583\ in-2 
I-O.577J ■ -"-U 


/+0.539N -|n-2 
1-0.606J ■ 


(±0.388) ■ 10-2 


mil) ■ 10-2 


1.11± 0.01 


800. - 1500. 


1071. 


2.597 ■ 10-2 


/+0.128^ . -|Q-2 
I-O.I39J 




(±0.098) • 10-2 


/+0.034\ -in-2 

1-0.031J ■ 


1.09± 0.01 


1500. - 2500. 


1901. 


7.208 ■ 10-3 


/+0.270^ . -|Q-3 
I-O.324J 


/+0.484\ 1 r|-3 
l-0.460J ■ 


(±0.263) ■ 10-3 


/+0.095N in-3 
1-0.080^ ■ 


1.07± 0.01 


2500. - 5000. 


3420. 


1.734 ■ 10-3 


/+0.051\ in-3 
1-0.067J ■ 


/+0.114N . -|Q-3 


(±0.070) ■ 10-3 


i'+0.022n 1^-3 
I-O.OI9J ■ 


1.05± 0.01 


5000. - 20000. 


8095. 


1.281 ■ 10-^ 


/+0.038\ in-4 
1-0.046^ ■ 


/+0.068N 10-4 
I-O.O57J ■ -"-U 


(±0.055) • 10-^ 


/+0.008N 1 n-4 
I-O.OI4J ■ -"-U 


1.04± 0.01 



Table 21: The QCD predictions for the differential dijet cross section as a function of Q^. Other details are as described 
in the caption to Table 12. 



Q2 range [GeV2] 


<Q2> 


R2+1 


Am- 


Aa. 


A cxp 
^Pdf 


A thco 
^pdf 


Chad ± AChad 


470. - 800. 


604. 


7.615 ■ 10-2 




/+0.443^ . 1Q-2 
I-O.475J 


(±0.122) ■ 10-2 


f+omo\ in-2 
I-0.017J ■ 


1.11± 0.01 


800. - 1500. 


1071. 


9.424 ■ 10-2 


/+0.539N -|r|-2 

Uo.sseJ ■ 


/+0.583N -|n-2 
1-0.586J ■ ^'^ 


(±0.122) ■ 10-2 


/+0.007\ in-2 

U0.025J ■ 


1.09± 0.01 


1500. - 2500. 


1901. 


1.148 ■ 10-1 


/+0.050N in-l 

K-ome) ' 


/+0.074N . -|Q-1 


(±0.015) ■ 10-1 


/+0.003N in-l 

1-0.004J ■ 


1.07± 0.01 


2500. - 5000. 


3420. 


1.380 ■ 10-1 


/+0.047N . in-l 
Uo.058J '-^ 


/+0.089N -|n-l 
I-O.O81J ■ 


(±0.022) ■ 10-1 


/+0.007^ . 1Q-1 


1.05± 0.01 


5000. - 20000. 


8095. 


1.829 ■ 10-1 


/+0.057N . 


1-0.096^ 


(±0.029) ■ 10-1 


1-0.012J 


1.04± 0.01 



Table 22: The QCD predictions for the dijet fraction as a function of Q^. Other details are as described in the caption 
to Table 12. 



□ 



(Q) [GeV] 


«s((g)) 






Aaf 




24.6 


0.1436 


0.0047 


+0.0032 
-0.0066 


+0.0045 
-0.0041 


+0.0105 
-0.0079 


32.7 


0.1396 


0.0048 


+0.0019 
-0.0036 


+0.0034 
-0.0034 


+0.0078 
-0.0059 


43.6 


0.1306 


0.0063 


+0.0041 
-0.0044 


+0.0023 
-0.0024 


+0.0054 
-0.0040 


58.5 


0.1276 


0.0075 


+0.0027 
-0.0041 


+0.0019 
-0.0016 


+0.0046 
-0.0033 


90.0 


0.1149 


0.0111 


+0.0056 
-0.0099 


+0.0009 
-0.0013 


+0.0037 
-0.0028 



Table 23: The ag values as determined from the QCD fit to the measured 
dijet fraction -R2+1 as a function of Q. For each bin in Q^, the mean value (Q), the 
extracted value of the strong coupling constant, ^^((Q)), the statistical uncertainty, 
Aaf*^*, the systematic uncertainty (not) associated with the energy scale of the jets, 
AaP (Aa^^^*), and the total theoretical uncertainty, AaJ*^, are given. 
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Figure 1: 



The measured differential dijet cross sections in NC DIS as functions of 
a) Zp^i, b) logiQ XBj, c) logj^Q ^ and d) dijet invariant mass Mjj. The inner error bars 
represent the statistical errors of the data. The outer error bars show the statistical 
errors and systematic uncertainties — except those associated with the uncertainty in 
the absolute energy scale of the jets — added in quadrature. For comparison, pure NLO 
QCD calculations (dashed lines) and NLO QCD calculations corrected for hadronisation 
effects (solid lines), obtained using the proton MBFIT PDFs and fiR = fip = Q, are 
shown. The relative differences of the measured differential cross sections over the NLO 
QCD predictions corrected for hadronisation effects are shown underneath each plot; 
the shaded band represents the uncertainty of the QCD calculation (see text). The 
hadronisation correction (Chad) together with its uncertainty are shown above each plot. 
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Figure 2: The measured differential dijet cross sections in NC DIS as functions of 
a) E^^^j^ , b) r/'g*'^, c) E^p"!^ and d) r/'g*'^. Other details are as described in the caption to 
Fig. % 
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Figure 3: a) The measured inclusive {datot/dQ'^) and dijet {da2+i/dQ'^) differential 
cross sections in NC DIS as a function of Q^. The hadronisation correction (Chad); shown 
above the figure, refers to the dijet cross section, b) The dijet fraction, i?2+i(Q^)) in NC 
DIS as a function of Q^. The light shaded band displays the uncertainty due to the 
absolute energy scale of the jets. For comparison, the QCD predictions using MBFIT 
proton PDFs determined assuming as{Mz) = 0.113 and 0.123 are also shown. The 
bands showing the theoretical uncertainty on the cross sections and dijet fraction do not 
include the uncertainty associated with as{Mz)- Other details are as described in the 
caption to Fig. ||. 
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Figure 4: The relative uncertainty on a) the inclusive, b) dijet differential cross 
sections and c) the dijet fraction, due to the statistical and systematic experimental 
uncertainties of each data set used in the determination of the MBFIT PDFs. The 
shaded (hatched) bands indicate the uncertainties obtained (not) taking into account the 
correlations among the PDFs parameters. 
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Figure 5: a) The as{Q) values determined from the QCD fit of the measured dijet 
fraction, R2+i{Q'^), as a function of Q. The inner error bars represent the statistical errors 
of the data. The outer error bars show the statistical errors and systematic uncertainties 
— except those associated with the uncertainty in the absolute energy scale of the jets — 
added in quadrature. The three curves indicate the renormalisation group predictions 
obtained from the PDG as{Mz) value and its associated uncertainty, b) The uncertainty 
on tts due to the absolute energy scale of the jets, Aaf^. c) The total theoretical 
uncertainty associated with the determination of a^, AaJ^. 
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